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ABSTRACT 

The continuity equations for neutral and charged particles in a partly ionized gas 
may possess periodic solutions. These oscillations are suggested as a possible clock 
mechanism for periodic astrophysical phenomena such as geomagnetic micropulsations, 
pulsars, and of long -period fluctuations in the charged particle density of the interstel- 
lar medium. A mathematical model of these oscillations is derived and applied to two 
cases of astrophysical interest; an isolated, virtually uniform region of partially ionized 
gas; and a slightly ionized stellar or planetary atmosphere, the characteristics of which 
vary with position. The theoretical analysis is supplemented by an experimental inves- 
tigation of the pulsed production and loss of electrons from a magnetically confined lab- 
oratory plasma which was subject to the continuity -equation oscillation. Some pulsar - 
like characteristics were observed and measured in this plasma. 
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POSSIBLE APPLICATIONS OF THE CONTINUITY -EQUATION PLASMA OSCILLATION 
TO PULSARS AND OTHER PERIODIC ASTROPHYSICAL PHENOMENA 

by J. Reece Roth 
Lewis Research Center 

SUMMARY 

The continuity equations for neutral and charged particles in a partly ionized gas may possess 
periodic solutions. These oscillations are suggested as a possible clock mechanism for periodic 
astrophysical phenomena such as geomagnetic micropulsations, pulsars, and long -period fluctua- 
tions in the charged particle density of the interstellar medium. Analytical and experimental 
investigations of this oscillation are discussed in light of their possible astrophysical relevance. 

A mathematical model is derived which exhibits such oscillations. The physical conditions that 
must occur in order that this model be a reasonable approximation are examined. 

The most restrictive of these conditions is that a population of energetic electrons, about 
100 electron volts in energy, must permeate the partially ionized gas in question; that the ioniza- 
tion mean free paths be larger than the minimum characteristic dimensions of the oscillating re- 
gion, and that the product of time -averaged electron and neutral number densities, n e n, be in- 
dependent of position within the coherently oscillating region. 

The particular model of these oscillations which is discussed is applied to two cases of as- 
trophysical interest. The first is an isolated, spatially uniform region of partially ionized gas, 
the characteristics of which are intended to be representative of interstellar H-I and H-H regions. 
The second case is one in which the electron and neutral number density may vary significantly 
with position, as would be the case in a planetary atmosphere or the upper reaches of a stellar 
atmosphere. It is shown that the physical conditions existing in these two cases may satisfy the 
conditions necessary to make the oscillation possible. 

Exact solutions to the nonlinear continuity equations were obtained on an analog computer. 

It was found that the electron number density waveform displayed the pulsar -like characteristic 
of sharp, narrow, periodic peaks between broad, flat minima. If a turbulent heating process or 
fluctuating electric fields were available to heat these pulses of electrons to relativistic energies, 
the pulses of electrons may well be accompanied by pulsed radio -frequency radiation. 

An experimental program was undertaken to study the pulsed emission of electrons observed 
from a laboratory plasma subject to the continuity -equation oscillation. Some suggestive charac- 
teristics have been observed and measured in this laboratory plasma, including the following: 

(1) long-term stability of the pulse frequency within limits set by apparatus limitations; (2) obser- 
vation of fine structure in the pulses of electron emission; and (3) the frequency of the bursts was 
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observed to vary over the range 87 < v < 10 hertz. Much lower frequencies are to be expected 
in an astrophysical context, as a consequence of the lower particle densities characteristic of 
such applications. 



INTRODUCTION 


Periodic pulses of radio -frequency (RF) emission have recently been observed ema- 
nating from sources in interstellar space (refs. 1 to 3). These pulses have been observed 
to have the following characteristics: 

(1) The pulses of RF emission have an extremely stable period, constant to one part 
in 10® or 10^ (refs. 1 and 2). 

(2) Several periodic sources have been observed, the pulse frequencies of which are 
about 1. 0 hertz. 

(3) This periodic RF emission has been received over frequencies from about 40 to 
600 megahertz (ref. 2). 

(4) The flux intensity appears to fall off with frequency according to a ® power 

law (ref. 2). 

(5) The RF radiation from these sources has been observed to be highly polarized 
(ref. 3). 



Pulsar CP0950 Pulsar CP1133 Pulsar CP1919 


Figure 1. - Typical consecutive pulse waveforms for three pulsars, 
taken from reference 6. Waveforms were observed at frequency of 
408 megahertz with bandwidth of 4 megahertz and time constant of 
1 millisecond. 
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Figure 2. - Output of two receivers at 408 megahertz which were connected to 
orthogonal dipoles, as reported in reference 3. Pulsar CP1133 was under 
observation. Note instances of strong polarization of received signal. 


(6) Periodic optical emission has not been observed from a pulsar, the presumed op- 
tical image of which has been studied (ref. 4). 

(7) The individual pulses of RF emission have a fine structure (refs. 5 to 7). 

(8) The leading edge of the RF pulse appears to have a rise time of 1 millisecond or 
less (refs. 5, 6, and 8). 

On figure 1 is shown a series of consecutive pulses, at a frequency of 408 megahertz, 
for three sources (taken from ref. 6, p. 327). On figure 2 is shown pulses at 408 mega- 
hertz from the source CP 1133, taken for two perpendicular components of radio fre- 
quency polarization (fig. 2 of ref. 3). 

When taken together, the above characteristics are quite puzzling, and cannot be 
readily explained in terms of a physical model the existence of which is known on the 
basis of independent observational evidence (refs. 2 and 9). It appeared to the author 
that the continuity -equation oscillation mechanism might be capable of providing a peri- 
odic, clock-like mechanism for the production of electrons. The observed RF pulses 
may result from the stochastic heating and synchrotron radiation of such pulses of elec- 
trons. 

Several entirely distinct physical models are capable of yielding a coupled pair of 
continuity equations that possess periodic solutions (ref. 10). A general model for as- 
trophysical applications is shown in figure 3. It is assumed that the oscillating region is 
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Neutral influx 



Figure 3. - Schematic of region of partially ionized gas sub- 
ject to continuity-equation oscillation. Such a region 
ionizes a small fraction of the neutrals streaming through 
it. Time delays associated with motion of neutrals across 
region give rise to time-dependent net efflux of charged 
particles from region, as discussed in appendixes B and C. 


populated by electrons sufficiently energetic to ionize the neutral gas within it. The 
ambient neutral gas can be conceptually divided into two categories: (1) Neutral gas that 
passes through the oscillating region without being ionized, shown on the top of figure 3, 
and (2) a small fraction of this neutral flux that interacts in the oscillating region during 
each cycle of oscillation. The electrons and ions produced in the oscillating region are 
either lost to space or recombine in the surroundings. Photoionization and three -body 
recombination processes are assumed either to be not present in this model, or, if pre- 
sent, to be in balance. 

In the model of figure 3 it is important to note that at any instant of time there is not 
a detailed balance between the oscillating flux of charged particles leaving the region of 
interest, and the virtually constant flux of neutrals entering it. The time delay associated 
with the transit of neutral particles across the oscillating region acts as an averaging 
mechanism that prevents a detailed balance at any point within the ionized gas, and makes 
possible a relaxation oscillation of the charged particle number density within it. This 
physical process will be elucidated in the subsequent discussion. 

The present author has analyzed a physical model similar to that of figure 3, which 
yields continuity equations identical to those to be derived below (refs. 10 and 11). The 
synchronous charged particle number density oscillations predicted by these analyses 
were observed in the author's laboratory under conditions that roughly simulate the as- 
trophysical case (refs. 12 and 13). The simultude with the astrophysical case became 
even more suggestive when it was observed that the discharge in question was emitting 
pulses of charged particles with the continuity -equation frequency. 
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THEORETICAL ANALYSIS 
Relevant Results of Prior Investigations 


Periodic solutions were obtained in reference 10 to the following pair of coupled, 
nonlinear differential equations 


n - C Q + Cjn + C 2 n g + Cgnn e + C^n 2 + 

(1) 

“e = A 0 + A l n + A 2 n e + A 3 nn e + C 4 n2 + C 5 n e 

(2) 




in the small amplitude limit, n(t) = h + 6n(t), where 5n(t) « n. The peak -to -peak am- 
plitude of n (t) may be large or small compared with the average value of n (t) . The 
coefficients A i and Cj of these equations are capable of representing the most com- 
monly encountered sources and sinks which appear in the neutral and charged particle 
continuity equations, with the exception of diffusion and three -body processes. 

These equations are written on the assumption that the Debye distance is smaller 
than the plasma dimensions, thus enabling one to eliminate either the ion or the electron 
continuity equation as an independent constraint on the oscillation. They also assume that 
the diffusion term is negligible, or may be represented by one of the coefficients or 

This assumption is justified in appendixes B and C. 

In reference 10 it was shown that the small -amplitude solution for n g (t) of the two 
equations (1) and (2) has a waveform given by Jacobian elliptic functions. These solutions 
are plotted on figures 4(a) to (c) for three particular cases. Figure 4(a) displays the 
sharp, isolated peaks separated by broad, flat minima that is characteristic of the Ja- 
cobian elliptic functions for large peak -to -peak amplitudes. 

In references 12 and 13 it was shown that the discharge under experimental investi- 
gation in these references could be adequately described by specializing equations (1) and 
(2) to the following form 


n = n n( crv) - n n(ov) _ 
e ' ne e ' ne 


n e = -nn e <ov> ne + nn e <ov) ne 


(3) 

(4) 


where n g and n are, respectively, the time and space averaged electron and neutral 
number density in the oscillating region. ( ov) is the ionization rate parameter for 
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electron -neutral ionization. In these references it was shown that the electron and neu- 
tral number densities are oscillating, with a frequency that is approximately 


v = 



(5) 


when the peak -to -peak amplitude of n g (t) is not comparable to its maximum value. The 
reader is referred to references 11 and 13 for a discussion of the case in which the peak- 
to-peak amplitude of n g (t) is comparable to its maximum value. In references 12 and 13, 
it is shown that the oscillations observed in an appropriate experimental situation obeyed 
the theory based on small -amplitude solutions to equations (3) and (4). 

These equations are not adequate if photoionization and three -body recombination 
processes are present in the region of interest, or if these processes are not in balance 
at all times. In some real cases, one might expect to find processes which demand con- 
sideration of the momentum or energy equations, or which involve time -dependent elec- 
tric or magnetic fields. The experimental results presented in reference 13 bear out the 
assumption that the electric and magnetic fields do not determine the frequency, ampli- 
tude, or waveform of the continuity -equation oscillation. 

In an astrophysical context, the second terms on the right-hand side of equations (3) 
and (4) represent, respectively, the instantaneous rate at which neutrals are lost, and 
charged particles are produced, by electron-neutral impact ionization. It is shown in 
appendixes B and C that a consideration of the divergence term of the continuity equations 
will yield the first terms appearing on the right-hand side of equations (3) and (4), pro- 
vided that the ionization mean free paths of electrons and neutrals are both larger than 
the minimum dimensions of the oscillating region, such that Fick's law is not valid. 


Requirements for Coherent, Periodic Oscillations 

In this section, the requirements will be discussed which must be met if this model 
of the oscillation process is to occur in an astrophysical context. The requirements which 
must be satisfied merely in order to make an analytical treatment possible will be dis- 
cussed in the next section. 

The basic physical process responsible for the oscillations under discussion is the 
time-dependent rate of ionization of neutrals by energetic electrons. The ionization-rate 
coefficient, (crv) ne , for this process is significant above an electron kinetic temperature 
of about 10 electron volts, and attains its maximum in the energy range from 80 to 400 


7 


mi Bill mm\ I 


electron volts. For neutral atoms of astrophysical interest, the ionization rate coeffi- 
cient varies by less than a factor of 2 from about 50 electron volts to 1 kiloelectron volts 
(refs. 14 and 15). In order for this ionization process to take place, there must exist a 
significant population of energetic electrons whose kinetic temperatures are above the 
ionization potential of the ambient neutral gas. 

An inspection of equation (5) suggests that a necessary condition for coherent oscil- 
lations due to this mechanism is that the density product n g n be constant in the region of 
coherent oscillation. The exact outcome of the continuity -equation oscillation for non- 
uniform n g n is not known. Presumably a synchronizing effect exists which can keep the 
oscillation coherent over minor variations from uniformity, since coherent oscillations 
were observed in the experiment described below when i^h varied by a factor of 10 over 
the discharge volume. If this product were to vary greatly with position, different regions 
may oscillate with different frequencies, and the oscillations in n g (t) would then be ob- 
served to be broad -band and incoherent. 

It is not difficult to see that the entire region can oscillate coherently if the electrons 
move across the characteristic dimension L of the oscillation region in a time less than 
the period of oscillation, T, 


However, this is not a necessary condition, as is discussed below in the section Synchro- 
nization of Pulsating Region. If the electrons have energy V g and velocity v g , we may 
use equation (5) to write equation (6) in the following form, 


— < - = — (7) 

V e 17 VV (av >ne 


For observed pulsar frequencies of about 1. 0 hertz, equation (7) suggests that the maxi- 
mum characteristic dimension L for 100 electron volts electrons must be less than 
6000 kilometers if this condition is to be satisfied. 

A necessary condition for the existence of synchronous oscillations is that the mean 
free time of the electrons in the oscillating region must be less than the period, if the 
peak-to-peak amplitude of the electron density is comparable to its maximum value. The 
electron mean free time is just 


lne= n{ov) 


ne 


( 8 ) 


8 



The requirement that this be less than the continuity equation period T is then ex- 
pressed by 


1 < 2xr 

n<<7V>ne < av >ne 

This may be rewritten in terms of the inequality 



This requirement is satisfied whenever the gas is partially ionized, such that n g 


( 9 ) 


( 10 ) 

«n. 


Requirements for Analytical Treatment 

A number of additional assumptions are desirable in order to facilitate an analytical 
derivation of the oscillations in question. The first such assumption is that the peak-to- 
peak amplitude of the fluctuations in neutral density is small by comparison with the 
average neutral density. This assumption also implies that n g « n Q , since there is an 
electron gained for each neutral lost. A second assumption that will be quite convenient 
during the course of the analysis below is that the electron or neutral velocity distribu- 
tion is monoenergetic or unidirectional. This should cause no difficulty in principle, 
since one can relate the assumed monoenergetic or unidirectional velocity to the most 
probable velocity of an isotropic, Maxwellian distribution. 

An attempt to apply the small amplitude approximation to continuity equations which 
include a diffusion term results in a nonlinear equation for which solutions are not known. 
One must at present restrict the analytical development to continuity equations without 
diffusion terms. This means that electrons and neutrals must not be subject to diffusion 
due to electron-neutral collisions, or that the mean free paths of electrons and neutrals 
for electron -neutral ionization must be comparable to or longer than the smallest charac- 
teristic dimension of the oscillating region. K v is the velocity of the neutral atoms, 
the mean free path of a neutral for an ionizing collision is 


X 


n 


v 

n e <ov) 


ne 


( 11 ) 
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The requirement that this be greater than the characteristic dimension L may be stated 


L < 


ri (ffv) 
e' ne 


( 12 ) 


The companion requirement that the mean free path of the electrons be large compared to 
the characteristic dimensions may be written 


L < 


n(crv) 


(13) 


ne 


where v g is the electron velocity. 

In the analysis in appendix B of a spherically symmetric H-I or H-II region, it is 
convenient to assume that the electron and neutral number densities are spatially uniform, 
except for the small oscillating spatial perturbations resulting from the periodic ioniza- 
tion process. The neglect of these spatial perturbations for pulsed sources of charged 
particles is shown to be justified by the excellent agreement of the theory based on this 
assumption with the experimental results reported in reference 13. 

In the study of planetary and stellar atmospheres in appendix C, it is convenient to 
assume that the constant product n g n arises from a neutral density that decreases ex- 
ponentially with distance above the surface, and an electron density that increases expo- 
nentially with distance above the surface, with the same scale height as the neutral atmos- 
phere. 


NATURE OF EXACT, OSCILLATORY SOLUTIONS 
Nondimensionalization of Equations 

It is shown in appendixes B and C that under the assumptions specified in the pre- 
vious section, the equations of continuity for a partially ionized gas of astrophysical in- 
terest may be approximated by a symmetric pair of equations, given by equations (3) 
and (4), 


" » n < S e - 
"e = - n e ( " 


n e K<7V> ne 

(3) 

-nXov> ne 

(4) 
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where n and n g are both functions of time, and may also be functions of position. n e 
and ii are the time and space averaged values of the electron and neutral number density. 

This pair of equations has long been known to possess oscillating solutions (refs. 16 
to 18). 

Exact closed-form solutions to these equations have not been obtained, because of 
their nonlinear nature. It is of some interest to obtain exact solutions numerically and 
compare the waveform geometry of n g (t) with, for example, that of the pulsar bursts 
shown in figures 1 and 2. If we apply the initial conditions 

h e (0) = 0 (14) 

n e (°) = n eo (15) 

to n g (t), equation (4) then implies that the initial condition on n(t) is 

n(0) = n (16) 


It should be noted that n g0 is an extremal value of n g (t), and may be either a maximum 
or minimum of the charged particle density. 

It is convenient for purposes of obtaining exact solutions to equations (3) and (4) to 
define a nondimensional neutral density, 

x 


= n(t) 


n 


(17) 


a nondimensional electron density 


an extremal fraction of ionization, 


n e (t) 

y s -^— 

n eo 


n 


(18) 


(19) 


and an ' ' amplitude index* * 


V 


e 



n 


eo 


( 20 ) 
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which is the ratio of the time averaged value of n g (t) to its extremal value. When n eQ 
is chosen as the maximum value of n g (t), then y ^ 1 . 0 , v e — 1 - 0 , and 6 is the maxi- 
mum fraction ionized. When n eQ is chosen as the minimum electron number density, 
y>1.0, 7 j g s 1. 0, and 5 is the minimum fraction ionized. It is also useful to define a 
nondimensional time, 


T =< n eo”> 1 / 2 < OT >ne t ( 21 ) 

which will make the frequencies and periods of order unity. 

Substituting equations (17) to (21) into equations (3) and (4), one obtains the continuity 
equations in dimensionless form, 


x = x 0? e - y)5 1//2 

( 22 ) 

y = -y(i - x)5 -1 / 2 

(23) 


where the dots now represent a derivative with respect to the dimensionless time t. 
This pair of equations can be solved numerically by a computer, or in closed form by 
using the small amplitude approximation. 


Analog Computer Solutions for n p (t) 

c 

Equations (22) and (23) were programmed on an analog computer, and the nature of 
the waveforms of y(t) studied as a function of the amplitude index 77 and of the maxi- 
mum fraction of ionization, 6 . 

The waveform of y(t) is illustrated in figure 5(a) through (k) for 8=0. 01. It is 
notable that when 77 g is near unity, the peak -to -peak amplitude of y(t) is small, and its 
waveform is sinusoidal. As 77 g approaches zero, however, the peak -to -peak amplitude 
of y(t) approaches unity, and the waveform takes on a characteristic appearance of sharp, 
narrow peaks, separated by broad, flat minama. These characteristics of the waveforms 
of y(t) for rf e ^ 0. 25 are particularly suggestive of the pulsar waveforms illustrated in 
figure 1. Also apparent on figure 5 is the fact that the frequency is a function of the am- 
plitude index rj e , and hence of the peak -to -peak amplitude of y(t). Such a dependence of 
frequency on amplitude is consistent with the nonlinear nature of equations (22) and (23). 

It was found that the peak-to-peak amplitude of the oscillation was not a function of 
the fraction ionized, 5. This observation is consistent with the predictions of the small- 
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Dimensionless electron number density, y{x) 






(k)T? e = 0.125. 
Figured - Concluded. 


amplitude solution obtained in references 10 and 11. A plot of the peak -to -peak ampli- 
tude of y(t), Ay is shown in figure 6 as a function of r) e on the curve labeled " exact. M 

In figure 7 is shown an exact analog computer solution for the waveform of y(t) as a 
function of 5 when ri e = 0. 25. This figure illustrates the small changes in frequency 
and geometry of the waveform that result from variation of the fraction ionized, 5, over 
the range 0. 04 ^ 6 ^ 1. 0. The curves for 6 = 0. 04 and 0. 10 are virtually coincident in 
this figure. It is evident that the peak -to -peak amplitude of y(t) is not a function of 6, 
and that the frequency increases slightly as 6 is varied from 1. 0 down to 0. 10. There 
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Figure 6 . - Peak-to-peak amplitude of oscillation of dimensionless electron number density as function 
of amplitude index r?. The exact results are taken from analog computer data and are not a function 
of maximum fraction ionized. Small-amplitude results fail below rj = 0.80 because of spurious 
limit introduced by small-amplitude approximation at 17 = 2/3. 


are also slight changes in the waveform of y(t) as 6 is varied over this range, with the 
peaks becoming more symmetrical about their maximum as 6 is reduced below 1.0. It 
is clear from figure 7 that there is no significant effect of 6 on either the frequency or 
the geometry of y(t) when 5 is below 0. 10. One may therefore safely ignore 6 as a 
significant independent variable when the maximum fraction of ionization of a gas is less 
than 10 percent. 

The nondimensionalization of time chosen in equation (21) allows the frequency of 
oscillation shown in figure 5(a) to (k) to be expressed in the form (ref. 11) 


(n eo fi ) 1//2 < ov > 


V = 


ne 


27T 


G(f? e , 8) 


(24) 


where G(tj , 5) accounts for the lowering of the frequency with rj shown on figure 5. The 
quantity G(tj , 6) is shown as a function of t] in figure 8 for the case 6=1.0 and 
6 = 0. 01. The latter curve is appropriate for all values of 6^0. 10, since G(?j e , 6) is 
virtually independent of 6 for 6^0. 10. 

A comparison is given in appendix D of these exact analog computer solutions with 
the small amplitude approximation of references 10 and 11. 
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Frequency reduction parameter, G(tj, 6 ) 



Figure 7. - Exact, analog computer waveforms of dimensionless electron number density for 
exact amplitude index of 0.25. Small effect of maximum fraction ionized 6 on shape of 
waveform and frequency of oscillation is shown. Note that peak-to-peak amplitude of os- 
cillation is not a function of 6 . 



Amplitude index, r} 

Figure 8 . - Frequency reduction parameter Gty, 6 ) as function of amplitude index tj and maximum fraction ion- 
ized 6 . Exact, analog-computer solutions show small effect of 6 on G. Below 6=0.10, effect is negligible. 
Value of G<r?) based on small-amplitude approximation fails below 17 = 0.80. 
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RELEVANT OBSERVATIONS IN A LABORATORY PLASMA 


Description of Apparatus and Discharge 

The superconducting magnet facility in which the experiments were conducted has 
been described previously in reference 19. The geometry of the discharge has been de- 
scribed in reference 20, and its properties and operation in references 21 and 22. An 
isometric cut-away drawing of the discharge and magnetic field coils is shown in figure 9. 
The plasma is a hot -ion, modified Penning discharge confined in a magnetic mirror 
geometry with a mirror ratio of 2. 5:1. The discharge is located in the center of a 0. 85 
meter diameter, 1. 7 meter long vacuum tank. The magnetic field lines on which the 
charged particles are confined intersect the tank walls about 50 centimeters outside the 
discharge. The magnetic mirror existing between the walls and the confining mirrors 
should act to prevent a backflow of charged particles from the walls to the discharge. The 
discharge was normally operated with a maximum magnetic field on the axis of 1. 0 tesla, 
and at pressures below lCT^ torr. At these pressures, the mean free path of electrons 
and neutrals for ionization was much larger than the discharge dimensions, a necessary 
condition for a meaningful simulation of astrophysical phenomena. 



Figure 9. - Isometric, cutaway drawing showing discharge in magnetic mirror, 
superconducting magnet facility. 
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Functional Dependence of Frequency on n e n 


The characteristics of this oscillation have been observed systematically and reported 
in detail in references 12 and 13. Data were taken which confirmed the functional depen- 
dence of frequency on the square root of the product n g n given by equation (24). A 
log -log plot of experimental data for deuterium gas, taken from reference 13, is shown 
in figure 10. The observed frequencies for this gas ranged from 340 to 50 500 hertz. 

The observed neutral density was based on the vacuum gage reading. Relative 
values of the electron density were measured by a technique based on analysis of the re- 
tarding potential curves from a gridded electrostatic analyzer, discussed in reference 22. 
The best -fitting straight line of slope 0. 5 is drawn through the data points in figure 10. 

The good agreement of this slope with the data, over five orders of magnitude in the 
product n n, confirms the square root dependence of v on the product n n given by 
equation (24). 



Figure 10. - Observed frequency of continuity equation oscillation as function of average electron and 
neutral number density product, in arbitrary units, for 247 data points taken with deuterium gas (ref. 13). 
The best-fitting straight line of slope 1/2 is drawn through data, confirming dependence of v on 
(n e n) 1/2 . 


Frequency Stability of Oscillation 

One of the outstanding characteristics of the pulsars is long-term frequency stabili- 
ty, with long-term variations being no more than one part in 10 6 or 10 7 during 6 months 
or more. 
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An experiment was made to assess the frequency stability of the laboratory plasma, 
and confirm that any drift observed could be accounted for by the observed drift in or 
measurement uncertainty of the electron and neutral number density of the plasma. In- 
dividual pulses of electrons reaching the collector plate of a multigrid electrostatic 
analyzer were counted for 30 seconds. A schematic diagram of the instrumentation is 
shown in figure 11. One such 30-second measurement of the frequency was taken at 
5 minute intervals over a 200 minute experimental run. Special precautions were taken 
to keep the experimental conditions as constant as the nature and limitations of the ap- 
paratus would permit. 

The result of this frequency stability run is shown in figure 12. During this 200 
minute run, the frequency varied by ±3 percent about its mean value. This variation is 
easily accounted for by the ±2. 7 percent variations that were measured in the neutral 
pressure, and the ±13 percent variations that were measured in the electron density. 
These apparent fluctuations in electron and neutral number density were probably largly 
due to imprecision of the measuring techniques used, amounting to ±2 percent for rela- 
tive measurements of the neutral density, and as much as ±20 percent for relative 
measurements of the electron density (ref. 22). The operating conditions could not be 
controlled sufficiently to test the ultimate frequency stabiliby of the oscillation. It is 
interesting to note, however, that the few percent variation of the 28 kilohertz frequency 
shown in figure 12 for 200 minutes implies that the frequency did not change by more 

Q 

than a few percent over about 3. 3X10 individual pulses. This would imply^a similar 
degree of frequency stability over about 10 years in a 1 second pulsar. 


Oscillating plasma 



Figure 11. - Schematic drawing of instrumentation used to study waveforms 
of pulsed electron output and light intensity. Photomultiplier integrated 
light intensity across diameter of discharge. 
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Figure 12. - Observed continuity-equation oscillation frequency as function of time, for a 200-minute experi- 
mental run. Variation of frequency was +3 percent of average value, which was within limits to which elec- 
tron and neutral number density could be controlled. 

Experimental Waveform of n e (t) 

The pulsar waveforms shown on figures 1 and 2 are of radio frequency emission, and 
the fine structure therein may be ascribed to the mechanism by which the RF emission 
is produced. Nonetheless, it is of some interest to measure the time -dependent electron 
density in the discharge, and see what similarities and differences exist with respect to 
the waveforms of the pulsars shown in figures 1 and 2. The waveform of the ion current 
flowing from discharge was measured by connecting the collector of the energy analyzer 
shown in figure 11 to ground through a 1 kS2 resistor. The voltage drop across this re- 
sistor was then displayed on an oscilloscope, and gives a measure of the ion current lost 
from the discharge. Since the Debye distance was small, the ion and electron currents 
flowing from the discharge should be equal and in phase. The oscilloscope traces of ion 
current should be virtually identical with the instantaneous electron and ion density within 
the discharge, except for a transit -time spreading of the ions with different velocities. 

A representative sample of the data showing ion pulses from the discharge is given 
in figure 13. The specific operating conditions under which each was taken is listed in 
table I. In each photograph, the time increases from right to left, and the zero level is 
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Figure 13. - Experimentally measured waveforms of pulsed electron flux from discharge 
subject to continuity-equation oscillation. Experimental conditions for each case are 
listed in table I. All runs shown were with deuterium gas, at a maximum magnetic 
field of B m ax= 1.0 telsa. Electron number density is proportional to anode current. 
Zero level is baseline, or lowest excursion, of trace between pulses. 
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(i) Run GJ 14. 


(k) | Unnumbered run. 



(I) Run GJ19. 



(n) Run GJ25. 




Figure 13. - Concluded. 


the baseline between pulses. The electron number density is proportional to the anode 
voltage shown in table I for a fixed neutral number density. In general, the pulses be- 
come broader and less regular as the neutral density was increased. The "tail" follow- 
ing the initial sharp pulse is not primarily due to a velocity selection effect, in which 
slower particles arrive after the faster. The transit time from the center of the dis- 
charge to the analyzer for an ion of average velocity was typically 1/10 the width of these 
tails. 

The light output was monitored with a type IP28 photomultiplier, the output of which 
was displayed on an oscilloscope. A representative sample of these records of light out- 
put is shown in figure 14. In these traces, the zero level is at the bottom of the graticule 
and time increases from right to left. 

The maxima in the light output was synchronous with the pulse of ions emitted from 
the discharge. The light and ion and electron efflux pulsations had the same frequency 
and were in phase. 
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TABLE I. - OPERATING CONDITIONS FOR WAVEFORMS OF 


n e (t) IN FIGURE 13 


Waveform 

Run 

Indicated 

tank 

pressure, 

torr 

Direct- 

current 

anode 

voltage, 

kV 

Average 

anode 

current, 

mA 

y-axis 
scale, 
mV per 
major 
division 

x-axis 
scale, 
|isec per 
major 
division 

13(a) 

GK2 

2.45X10" 6 

10 

3 

1. 

0 

200 

13(b) 

GK3 

2.45 

15 

4.5 



200 

13(c) 

GJ2 

4.75 

10 

4 



50 

13(d) 

GJ6 

8.1 

5 

2.5 



50 

13(e) 

GJ8 

8.1 

15 

8.5 



2 

1 

13(f) 

GJ10 

8,1 

25 

14 





13(g) 

GJ11 

1.2X10" 5 

5 

3. 5 





13(h) 

GJ13 

1.2 

15 

11.5 





13(1) 

GJ14 

1.2 

20 

15 



10 

13(j) 

GJ18 

1.65 

15 

17 

2. 

0 

20 

13(k) 

— 

2.2 

10 

14 

2. 

0 

20 

13(1) 

GJ19 

1.65 

20 

21 

2. 

0 

10 

13(m) 

GJ22 

2.6 

10 

24 

5. 

0 

10 

13(n) 

GJ25 

2.6 

25 

46 

5.0 

10 


TABLE n. - OPERATING CONDITIONS FOR WAVEFORMS OF LIGHT 
INTENSITY IN FIGURE 14 


Waveform 

Rim 

Indicated 

tank 

pressure, 

torr 

Direct- 

current 

anode 

voltage, 

kV 

Average 

anode 

current, 

mA 

y-axis 
scale, 
mV per 
major 
division 

x- axis 
scale, 
fi sec per 
major 
division 

Gas 

14(a) 

Gi9 

2.6X10" 6 

4 

1.3 

10 

100 

He 

14(b) 

Gi8 

2.65 

5 

1.6 

10 

100 



14(c) 

Gi6 

2. 8 

10 

3.3 

20 





14(d) 

Gi2 

3. 1 

20 

7.1 

100 





14(e) 

CR10 

7.4 

5 

3 

200 



D 

2 

14(f) 

CR12 

7.4 

8 

5 

200 



d. 

2 





(e) Run CR 10. + 1 (f)RunCR12. 


Figure 14. - Light output of discharge for experimental conditions given in table II. 
Zero level is the bottom of graticule. 


QUANTITATIVE ASPECTS OF ASTROPHYSICAL APPLICATIONS 
Application to Interstellar H-I and H-II Regions 

One possible astrophysical application of equations (3) and (4) is the behavior of 
interstellar H-I and H-II regions, in which there exists a population of electrons suffi- 
ciently energetic to ionize any neutral hydrogen that may be present. The periodic fluc- 
tuations in electron density that would arise from this mechanism may be observable 
indirectly as long-term oscillation of point radio sources, or periodic fluctuations in the 
angle of rotation of polarized radio waves from a background radio source. If the period- 
ically fluctuating electron population is increased greatly in energy by interstellar elec- 
tric fields or plasma turbulence, then periodically fluctuating radio frequency emission 
may result from bremsstrahlung or cyclotron radiation mechanisms. 

The expected period of such oscillations in the interstellar medium is approximately 
given by equation (27). 
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27 T 


(25) 



ne 


The reaction rate parameter (crv) ng for atomic hydrogen has a maximum value at 
100 electron volts electron energy given by 


(ov) 


ne max 


3xl0 -8 cm^/sec 


(26) 


(ref. 15, fig. 25). The accepted maximum value of the interstellar neutral hydrogen 
density is n«l atom per cubic centimeter, and a charged particle number density of 
n g x 0. 5 electron per cubic centimeter (ref. 23, p. 254). These values, when substituted 
into equation (25) with the reaction rate parameter in expression (26), yield a period of 


T 


~ 277 V2 
3X10" 8 


« 10 years 


(27) 


If the electron population were cooler or substantially hotter than 100 electron volts, 

(crv) ne would be smaller than that given by equation (26), and the period would be larger 
than indicated by equation (27). These long -period fluctuations in the density of the inter- 
stellar electron number density may give rise to scintillations in the strength and polari- 
zation of radiation from radio stars. Such scintillations have not been observed, but this 
may be because radio -astronomical data have not been taken for a sufficiently long time. 

If the oscillating region had a higher electron and neutral density, the oscillation 
periods may be on the order of months. Oscillations due to this continuity -equation mech- 
anism may be responsible for the quasi -periodic fluctuations, with characteristic times of 
the order of months, that have been observed in the X-ray emission of certain astronomi- 
cal objects (refs. 24 and 25). 


Application to Planetary Atmospheres 

From the examples studied to date, it appears that the upper reaches of planetary 
atmospheres contain both neutral and charged particles. The charged particle density is 
much enhanced if the planet in question has a dipole magnetic field strong enough to trap 
such particles. 
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It is shown in appendix C that in-phase oscillations with a single frequency are pro- 
bable if the product n e n is constant over a significant region in the planetary atmosphere. 
A more usual situation is one in which n e n varies with height. Each height probably 
would oscillate in electron number density with a characteristic frequency of its own, and 
the net result would then be a broad spectrum of electron number density oscillation fre- 
quencies. 

The presence of a planetary magnetic field greatly enhances the likelihood of the os- 
cillations of n g (t) having observable effects. The periodic loading of the geomagnetic 
field lines with electrons may, for example, result in observable low-frequency periodic 
fluctuations in the strength of the magnetic field at the planetary surface. Such low- 
frequency oscillations - called "geomagnetic micropulsations" - are observed on the sur- 
face of the earth (refs. 26 to 28). Their origin is obscure, and they may be caused by 
continuity -equation oscillations in the electron number density of the magnetosphere. 

If the electron number density is oscillating with a single frequency over a signifi- 
cant vertical distance, if a strong geomagnetic field is present, and if the electrons are 
raised quickly to higher energies by turbulent heating or some other process, one might 
then expect to find radio -frequency emission from the trapped electrons, modulated by the 
continuity -equation frequency. It is very doubtful, however, whether such planetary 
"pulsars" have enough energy available for RF emission to be detected outside their own 
solar system. 

It is of interest to run a quantitative check and see whether the electron and neutral 
number densities in the earth's exosphere are consistent with the geomagnetic micropul- 
sation frequencies actually observed. It is known that there is no extensive region in the 
earth' s magnetosphere in which the product n g n is constant; it decreases monotonically 
with radius above an altitude of 100 kilometers. One would then expect a continuous dis- 
tribution (spectrum) of frequencies. Such a continuous distribution is characteristic of 
the geomagnetic micropulsations. The period of the geomagnetic micropulsations, if 
given by equation (25), is inversely proportional to (n n)^^. This implies that the period 
of the pulsations ought to be relatively long during the quiet times of solar sunspot mini- 
mum, when n is low. Near sunspot maximum, when n is enhanced by large factors 
over the values characteristic of sunspot minimum, one would expect the geomagnetic 
micropulsation period to be shorter. Such a trend of shorter period micropulsations 
with increasing solar activity (and hence indirectly with increasing n e in the magneto- 
sphere) has been observed (ref. 26). 

In reference 28 it is reported that geomagnetic micropulsations were observed simul- 
taneously at conjugate points in the geomagnetic field; observations of in-phase oscilla- 
tion over distances as great as 1000 kilometers are quite common (ref. 28). This sug- 
gests that, if the continuity -equation oscillation is responsible for the oscillations, they 
must originate at a sufficiently high altitude such that the electron mean free path is about 
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10 kilometers. On figure 15 is shown a log-log graph of the electron mean free path, 

and electron and neutral number densities as a function of height above the earth’s sur- 

3 ^ 

face. The electron mean free path is 10 kilometers at an altitude of about 500 kilome- 
ters. At this altitude, the electron number density is approximately 4xl0 5 per cubic 
centimeter, and the neutral density is approximately 5X10 per cubic centimeter. The 
dominant species at an altitude of 500 kilometers is atomic oxygen (ref. 29), whose 
maximum ionization rate parameter is given by Lotz (ref. 15, fig. 32) as (ctv) « 

^ ^7 no y 

10 cubic centimeter per second. Substituting these values into equation (32) we obtain 
a continuity -equation period of 

T « — « 14 sec (28) 

10“ 7 (4X10 5 X5X10 7 ) 1//2 

This frequency is quite close to the maximum spectral amplitude of about 0. 05 hertz 
(period 20 seconds) actually observed in the geomagnetic micropulsations, a spectrum of 
which is shown in figure 16 (ref. 27). 


Application to Stellar Atmospheres 

The preceding discussion of continuity -equation oscillations in planetary atmospheres 
is relevant in its qualitative features to stellar atmospheres, with two important diff er - 
ences: The first difference is that near a stellar surface, large amounts of energy are 
available for conversion to radio -frequency power. Enough energy is available so that if 
a conversion process is active over an entire stellar surface, even a very inefficient 
process can radiate enough power to be detected at interstellar distances. The second 
difference is that if pulsars are to be explained in terms of the periodic generation of elec- 
trons by the continuity -equation oscillation mechanism, then there must exist a large re- 
gion in the stellar atmosphere over which the product n g n is constant. 

The decrease of neutral density with height above a stellar surface may be approxi- 
mated by the isothermal law of atmospheres, 


n(Z) = n^e 


-Z/H 0 


( 29 ) 


The requirement of in-phase oscillations of a single frequency then requires that 
n(Z)n (Z) be constant, or that n (Z) have the functional form 

6 c 
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( 30 ) 


+Z/H 

n (Z) ~ n e 

e K ’ e, max 

where H e =H Q and the electron number density increases with distance above the stellar 
surface. This increasing n g could be the base of a stellar magnetosphere, in which the 
electron population increased exponentially as one moved from the neutral atmosphere in- 
to the base of the stellar "Van Allen Belts. " Such a stellar magnetosphere would imply 
a strong magnetic field, which would promote radio -frequency emission from energetic 
electrons. 

Substituting equations (29) and (30) into equation (25), we have an expression for the 
period of a pulsar, 


T = 


2ir 


2. 1x10 


8 


(fi e 6) 1/2 


<CTV> 


ne 


(n e n) 1/2 


sec 


(31) 


where the maximum ionization rate parameter for hydrogen, in equation (26), has been 
used. For the pulsar CP1919, with a period of 1. 337 seconds, equation (31) implies a 
product of electron and neutral number density. 


n e n 


2. 46xl0 16 /cm 6 


(32) 


This appears to be a reasonable number density product to be found in a stellar atmos- 
phere (ref. 23, p. 161). 

If the continuity -equation oscillation is confirmed as governing the frequency of pul- 
sars, a measurement of their period would then yield information about the product of 
electron and neutral number densities in a large region of the stellar atmosphere. 

It may be objected that the great stability of the pulsar frequency is inconsistent with 
the continuity -equation oscillation, since the electron and neutral number densities may 
be expected to fluctuate significantly in a typical stellar atmosphere. While this is true 
locall y, the expression for n(Z) given in equation (29) involves only quantities whose 
values averaged over an entire stellar surface (n^, H Q ) would be expected to vary only on 
an astronomical time scale. Temporal constancy of the electron number density is under- 
standable on the hypothesis of a very strong dipole magnetic field, which may average the 
stellar magnetospheric production and loss processes over the entire stellar surface, and 
over trapping times much larger than a pulsar period. Such an averaging process should 
be effective if the star is not subject to the equivalent of our 22 -year solar sunspot cycle, 
or if the number of charged particles stored in the stellar Van Allen Belts is much larger 
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than the number gained or lost over the length of time during which the oscillation period 
is observed to be constant. 


DISCUSSION 

Synchronization of Pulsating Region 

The pulsar waveforms shown on figures 1 and 2 have rise times as short as one mil- 
lisecond. It is usually argued that this rise time, when multiplied by the speed of light, 
places an upper bound on the size of the oscillating region. This limit results from the 
implicit assumption that a triggering signal must propagate throughout the volume of the 
oscillating region during the fast rise portion of the waveform. 

It is a highly significant fact that, if the continuity -equation oscillation is responsi- 
ble, such a triggering signal is not required for the production of synchronized pulses, 
and the size of the pulsating object is not limited to the product of the rise time and the 
speed of light. To see how this may be so, consider a "thought experiment" performed 
with two identical clocks, synchronized to one part in 10 . Suppose that these two clocks 
emit synchronized pulses every second, with an initial rise time of 10"° seconds. If 
these two clocks are placed on the opposite limbs of a star separated by a distance of 
(say) 1 light-second, an observer in an external stellar system will see the situation 
shown in figure 17. The observer will note the 1 microsecond rise time of the signals, 
and, being unaware of the situation, will incorrectly conclude that the pulsating object is 
no more than 300 meters in diameter, whereas in fact the clocks are 1 light second, or 
300 000 kilometers apart. 

It is important to realize that, if the ionization rate parameters (crv) are the same, 
two spatially separated regions of partially ionized gas will behave like two independent 
clocks whose frequency of oscillation is determined only by the product n g n. An exami- 
nation of the laboratory data presented in figure 13 shows that the rise time of the electron 
number density can be as small as 1 microsecond in this laboratory plasma, so the 
continuity -equation oscillation mechanism is capable of producing pulses of electrons 
with very short initial rise times. 


Distant 

observer 



Figure 17. - Observer looking at two synchronized clocks separated by 1 light second. 
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Figure 18. - Schematic of two adja- 
cent volumes in region subject to 
continuity-equation oscillation. 

If diaphragm is removed at time 
t - 0, mixing will occur, and 
products of electron and neutral 
number densities will become the 
same. Thus, both regions will 
oscillate with the same frequency. 
Phases of oscillation can become 
synchronized by mechanism dis- 
cussed in text. 


One can understand how adjacent regions can assume the same frequency and phase 
by considering the schematic arrangement in figure 18. Suppose that initially the two 
adjacent regions are isolated and have different values of n g n. If the diaphragm is re- 
moved at time t = 0 any transport process such as collisions, turbulence, ionospheric 
winds, etc. will cause a mixing of the electron and neutral number densities in adjacent 
regions. This will lend to both regions having the same values of h e h and the same 
oscillation frequency. 

It is also understandable that the two regions should have the same phase, as well 
as the same frequency. If the right-hand region pulses slightly in advance of the left- 
hand region, the particles which are lost from the right-hand region will pass over into 
the left-hand region and trigger it to pulse before it would have done so without an outside 
disturbance. The particles which move from one region to surrounding regions therefore 
can act as a synchronizing and phase -locking mechanism. One can use an obvious exten- 
sion of the above argument to show that after an initial period of incoherence and adjust- 
ment, a region of very large dimensions may be synchronized, if the product n e n is 
constant throughout the region. 


Possible Experimental Checks 

Any observations of a pulsar may reveal fine structure in the RF radiation which is 
consistent with many small regions within the source radiating independently, but with 
the same modulation frequency. These regions would have been excited by an expanding 
spherical shell of electrons from the oscillating layer. The small phase differences ob- 
served between individual peaks in the fine structure can be attributed to time delays at 
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the source arising from RF emission at different places and times near the pulsar, or to 
pulses of electrons originating at different locations within the oscillating region. Fine 
structure therefore need not be a characteristic of solutions to the continuity equations, 
but is consistent with probable secondary processes which produce the radiation. The 
existence of fine structure, recently reported by Craft, et al. (ref. 7), is consistent with 
this model. According to Craft, et al. (ref. 7), if the pulses shown on figures 1 and 2 
are examined under high temporal resolution, they resolve themselves into hundreds or 
thousands of individual sharply peaked pulses. 

When the periods are comparable or greater than the electron transit time across 
the oscillating region (as is the case with geomagnetic micropulsations), one would then 
expect the source to behave as a single oscillating region, with no fine structure to the 
pulses. 

An experimental check can be provided by an observation of light output. The reac- 
tion rate parameter (c rv) for excitation of atomic hydrogen is of the same order of mag- 

nitude as the ionization rate parameter (crv) ne , and one would therefore expect oscilla- 
tions in the electron number density to be accompanied by in-phase oscillations of the 
light. This was, of course, observed in the laboratory experiment discussed above. The 
peak -to -peak fluctuation in visible light will be quite small by comparison with the radiant 
intensity of a stellar surface, unless the star in question is extremely efficient in convert- 
ing its total energy output into energetic charged particles. One would, therefore, not 
expect to be able to observe fluctuations in light output when the oscillating source is in a 
stellar atmosphere. If such oscillations were observed, they should be of the same fre- 
quency and phase as the fluctuations of electron number density. If observations are 
confirmed of fluctuatuons in light output with a frequency different from the pulsar fre- 
quency, this would tend to rule out the continuity -equation oscillation mechanism, unless 
complicated additional physical mechanisms are invoked. 

When there is no bright background source behind the pulsating region, it may then 
be possible to detect long -period fluctuations in the light output of the pulsating region 
which would be in phase with the fluctuations in electron number density. Such an obser- 
vation does, of course, presuppose coherent oscillations and an adequate optical depth of 
gas. 

The pulsars observed thus far all seem to posess characteristic sharp, narrow 
peaks separated by broad, flat minima, whose minimum values are at or below the radio - 
frequency noise level. If the pulses are due to the continuity -equation oscillation, then 
this waveform suggests that r\ of equation (20) has either values n » 1.0, or 
77 e « 1. 0 (ref. 11). One might expect to find, sooner or later, a pulsar for which 
?j e « 1. 0. The radio -frequency waveform in this case would be expected to have a sinu- 
soidal modulation, and a constant, nonzero minimum level of radio -frequency emission, 
well above the background noise level. 
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K independent means are available to measure the electron and/or neutral number 
densities, one should observe a correlation between changes in the oscillation period and 
changes in the products n g n. If the period is constant (as it appears to be in pulsars), 
then this product should be observed to remain constant. 


CONCLUSIONS 

The continuity -equation oscillation provides a simple, straightforward physical 
mechanism by which electrons can be produced in periodic pulses by a partially ionized 
plasma. 

Arguments have been put forward to show that this mechanism is quantitatively and 
qualitatively consistent with such periodic astrophysical phenomena as pulsars, geomag- 
netic micropulsations, and periodic variations in the electron number density of interstel- 
lar H-I or H-II regions. 

The continuity equations have been derived for a grossly uniform spherical mass of 
gas, and for a planetary or stellar atmosphere in which the neutral and charged particle 
number density vary with height above the surface. It has been shown that if the ioniza- 
tion mean free paths are longer than the smallest dimension of the region, an oscillation 
in electron number density will occur which arises from transit time delays associated 
with the passage of neutrals across the ionized gas. These transit time delays are signi- 
ficant in that they prevent a detailed balance between the production and loss terms of the 
continuity equations. 

If the ionization mean free paths are long compared to the characteristic dimensions 
of the ionized gas, and if energetic electrons are available to cause ionization, then peri- 
odic oscillations in electron number density are a natural outcome of the operation of the 
continuity equations. Such oscillations have been observed over more than five orders of 
magnitude in the product n e ii in the laboratory, and frequencies as low as 87 hertz have 
been experimentally observed. Much lower frequencies would be expected in the lower 
electron and neutral number densities characteristic of astrophysical conditions. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio October 20, 1968, 

129-02-08-06-22. 
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APPENDIX A 


SYMBOLS 

[All quantities are in SI units unless noted otherwise] 

coefficients of terms in charged particle continuity equation 
coefficients of terms in neutral particle continuity equation 
G(? 7 , 6) frequency reduction parameter 

H e scale height of electron number density, assumed to increase exponentially 

with height 

H q scale height of neutral, isothermal atmosphere 

L characteristic dimension of oscillating region 

n neutral density 

n time and space averaged value of neutral particle density 

n time averaged neutral number density, in region where n is virtually constant 

in space 

n g electron density 

n 0 time and space averaged electron number density 

n time averaged electron number density, in region where n is virtually con- 

stant in space 

n „„ time -dependent electron number density at top of oscillating region in a plane- 
g ? rn ctx 

tary or stellar atmosphere 

n initial extremal value of n ft) at t = 0 

eo e ' 

n Q surrounding neutral gas density 

n.. neutral number density at Z 

nj time averaged neutral number density at Z m ^ n 

R radius of spherical oscillating region 

r- radial coordinate 

T period of oscillation, sec 

t time 

t d transit time of downward -moving neutral particle 
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mean free time between electron-neutral collisions 

transit time of neutral across oscillating region 

transit time of electrons across characteristic discharge dimension 

transit time across radius R of spherical region 

transit time of upward-moving neutral particle 

electron energy, eV 

neutral velocity 

electron velocity 

dimensionless neutral density, defined by eq. (17) 
dimensionless electron density, defined by eq. (18) 
height above surface of a star or planet 

upper boundary of oscillating region in a planetary or stellar atmosphere 
lower boundary of oscillating region in a planetary or stellar atmosphere 
argument of electron number density variation, defined by eq. (E3) 
phase angle resulting from transit time across discharge, defined by eq. (E4) 
extremal fraction of ionization 

small time -dependent perturbation of neutral density 

space and time averaged attenuation coefficient of neutrals, averaged over 
prior path in discharge, defined by eq. (B4) 

time averaged attenuation coefficient of downward -moving particles 

time averaged attenuation coefficient of upward-moving particles 

electron amplitude index used in exact numerical solutions of eqs. (22) and (23), 
defined by eq. (20) 

amplitude index used in small -amplitude solutions 

mean free path of a neutral for an ionizing collision 

frequency, Hz 

net neutral flux at a point 

net flux of electrons through a point 

total neutral particle flux through a point 

flux of upward -moving neutral particles 
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i// dimensionless attenuation coefficient, in spherical region, defined by eq. (B5) 

dimensionless attenuation coefficient of downward -moving particles, given by 
eq. (C9) 

j// u dimensionless attenuation coefficient of upward-moving particles, defined by 

eq. (C7) 

p path length of neutral in spherical oscillating region 

<CTv) ne ionization rate parameter, for ionization of neutrals by electron impact 

r dimensionless time, defined by eq. (21) 

9 polar angle in spherical coordinates 
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APPENDIX B 


DERIVATION OF CONTINUITY EQUATIONS IN A SPHERICALLY 
SYMMETRIC ISOLATED REGION OF PARTIALLY IONIZED GAS 

One possible application of the continuity-equation oscillation mechanism is to par- 
tially ionized interstellar H-I or H-II regions. The spherically symmetric mass of gas 
shown in figure 19 will be adopted as a model of such a region. Both the electron number 
density and neutral density are assumed uniform over the region, except for the very 
small spatial and temporal perturbations caused by the ionization process itself. The 
electrons may be confined to the region of interest, or may also exist in the region out- 
side the discharge. The electrons could, for example, be confined to the region of in- 
terest by a local magnetic field perturbation in the form of a magnetic bottle. For this 
analysis, it is assumed that the neutral atoms are isotropic in velocity space, and mono- 
energetic. 

Consider the net flux of particles at a position r from the center of a spherical 
region of radius R. The differential flux of neutrals reaching position r has been at- 
tenuated by an amount proportional to the path length p in the region. If the ionization 
mean free path of a neutral is much larger than the characteristic dimension of the re- 
gion, L « 2R, the radial flux through a unit area at the point r with a normal pointing 
along r is 


d< ?n (e) = \ n o ve " eP 


cos 9 sin 9 d 9 


(Bl) 


Incoming 



Figure 19. - Geometry of spherical oscillating region, assumed to be 
surrounded by neutral gas having an isotropic, monoenergetic 
velocity distribution. Only neutral particles whose velocity vectors 
have a polar angle 0 will passthrough point r. 
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where one has used spherical symmetry to integrate over the azimuthal angle. The re- 
gion is assumed surrounded by neutral gas of density n Q . This net flux is zero for an 
isotropic velocity distribution. The instantaneous value of the absorbtion parameter e 
is the inverse of the electron -neutral mean free path, 


e 



V 


(B2) 


This must be averaged over all previous times during which a particle passing through 
the point r could have been in the discharge, and hence subject to ionization. The max- 
imum length of time a neutral particle can be in the discharge before reaching the point 
r is 


t 


o 


= P 
v 


(B3) 


The average absorbtion coefficient is found by integration over all times from t-t Q to t, 


e = 



“e^X^ne 


df = 


v 


( <rv) 

2* / n (tf)dt' 

0 


(B4) 


since ( CTV ) ne and v are not functions of time by hypothesis. It is helpful to define the 
attenuation coefficient 


xl/ =ep = 


<av) 


ne 


f 

J t-t 


n e (t')dt' 


(B 5) 


The differential flux of equation (Bl) may then be written 

d cp (0) = - n ve~^ cos 0 sin 0 d0 (B6) 

n 2 ° 

It should be noted that \js is a function of both r and 0, through its dependence on p 
in the lower limit t . The net radial flux at position r is given by integrating over all 
possible polar angles. 
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</> n ( r ) = - n n v 



sin 9 cos 9 d0 = n v , 
o d 


(B7) 


where v^ is the macroscopic drift velocity of neutrals into an elemental volume around 
the point r. 

The neutral particle continuity equation can be written 

n= -V ■ cp n - nn e <ov) ne (B8) 


The divergence term in this equation can be written in spherical coordinates, by using 
equation (B7) 



The derivative within the integral sign in equation (B9) may be obtained by using equa- 
tions (B3) and (B5) which yield 


d^_ ^^pX^ne dp 
dr v dr 


(BIO) 


By applying the law of cosines to figure 19 one can obtain a relation between p, r, 
and cos 9, as follows: 


2 „2 2 

p = R -r + 2rp cos 9 

differentiating this with respect to r, one obtains 

dp p cos 9 - r 
dr p - r cos 9 


(Bll) 


(B12) 
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Substituting equations (BIO) and (B12) into equation (B9) one obtains for the divergence 
term 


n v r v 


e” ^ cos 0 sin 0 d0 


n ( ov) 
o' ne 


/ 7T 

n (t-t ) (P cos 9 . ~ .. r) e -1 / 7 cos ^ gj[ n $ (B13) 

e 0 (p - r cos 0) 


To simplify the analytical and physical aspects of the problem, we specialize to the 
position r = 0 at which if/ is no longer a function of 8, and equation (B3) becomes 

t = (B14) 

0 v 

By evaluating the second term in equation (B13) at r = 0, and applying l ' Hospital’s 
rule to the first term, one obtains the value of the divergence term at the center of the 
spherical mass of gas, 


-V ■ 


<p = n n (t-t )(<tv) e 

^n o e' o ' ' ' ne 


-if/ 


(B15) 


It is now possible to use a calculation of the total flux at r = 0 to calculate the time- 
dependent flux there. The total flux may be written by analogy with equation (B7) 


cp T ( 0, t) = n(0, t)v = 


n o ve 


-if/ 



sin 8 cos 9 dd 


(B16) 


from which it follows that the time -dependent neutral density at the center of the oscillat- 
ing region may be written 


n(0,t) = n o e"^ (B17) 

Substibuting equations (B15) and (B17) into equation (B8) for the neutral particle continuity 
equation yields 


(B18) 
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n=n(av>ne[ n e ( t - t o) " n e W ] 



It is clear that n is not zero because of the time lag represented by the retarded time 
t = R/v. In practice, the velocity of the neutrals will not have a constant value v, but 
will have a Maxwellian spread of energies, all the way down to zero velocity. Since some 
particles travel with extremely small velocities, these will average the attenuation over 
one or more periods of oscillation. We can then replace the retarded electron number 
density n e (t -t Q ) by the time averaged number density n e to yield 

n= n(ow> ne (n e -n e ) (B19) 

The validity and limits for which one can replace n e (t-t Q ) with n e are justified by exper- 
imental results discussed in appendix E. 

One can obtain the continuity equation for the electrons by an argument virtually iden- 
tical with that above. The electron continuity equation may be written 

n = -V • c? + nn (crv)„ Q (B20) 

e ^ne e ne ' ' 

The divergence term can be obtained by replacing n (t) — n(t), n(t) — n (t), etc. , and 

r*J C C 

realizing that the argument of the exponential e e P in equation (B16) is positive, since 
the ionization process always results in an excess of electrons in the oscillating region. 

By following a line of argument like that from equation (Bl) to equation (B19) one can show 
that the electron continuity equation can be written 


n e = n e <ov) ne [n(t) - n(t-t Q )] (B21) 

One can replace the retarded value of n(t-t Q ) (with t Q now defined in terms of the 
electron velocity v g rather than the neutral velocity v) by appeal to the same arguments 
cited for the previous case, 


n (t-t 0 ) « n 

so that equation (B21) assumes the symmetrical form 

n e =n e <av> ne (n -^ 


(B22) 


(B23) 
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APPENDIX C 


DERIVATION OF CONTINUITY EQUATIONS IN AN OSCILLATING 
REGION IN A STELLAR OR PLANETARY ATMOSPHERE 

In a planetary or stellar atmosphere, the electron and neutral density will be depend- 
ent upon the distance above the surface. It is, therefore, necessary to generalize the 
derivation in appendix B to the case for which n and n g are functions of position as 
well as time. A model of such an atmosphere is shown on figure 20. The planetary or 
stellar surface is located at Z = 0. The oscillating region is located between a lower 
boundary at Z m ^ n and an upper boundary at Z max . The neutral density in an isothermal 
atmosphere may be written 


n = n^(t) exp 


< Z - Z min>' 


H. 


o J 


(Cl) 


where H Q is the scale height of the atmosphere. 

It is assumed that the neutral particles viiich move between Z min < Z < z max 
are monoenergetic with the velocity v, and that half of these neutrals are directed upward, 
and half in the downward direction. It is assumed that the mean free path of both the 
electrons and neutral particles in this region are much longer than the vertical dimen- 
sion of the region itself, so that the probability of a neutral being ionized is much less 
than unity for a single traversal of the oscillating region. 



Figure 20. - Geometry of oscillating 
layer in planetary or stellar at- 
mosphere, confined between heights 
Z m in and Z max above surface. 
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A downward -moving neutral particle at the height Z will have been in the oscillating 
region for a transit time t^ given by 



Z -Z 
max 


and an upward moving particle for a transit time t , 


(C 2) 


t 


u 


Z-Z 


min 


(C3) 


The neutral flux will be attenuated exponentially as it moves through the oscillating region. 
The upward flux may be written 


<P U = “ n i v ex P 
2 


z-z 


mm 


H 


o J 




(C4) 


where the upward attenuation factor for neutrals must be written in terms of an average 
over the time and space varying conditions along the neutral path length. The instanta- 
neous absorbtion coefficient is, by analogy with equation (B2) 


e u = 


"e^t’K^ne 


(C 5) 


If the neutral particle is located at position Z at time t, its position Z^ at the 
prior time t' may be written 


Z^= Z - v(t-t’) 


(C 6) 


When averaged over the time which the neutral particle has been in the oscillating region, 
the attenuation coefficient \[/ u is 


xjs = (Z-Z . ) € = (Z " Zmin) 
v min' u 


t 


/ 

J t-t 

L l u 


t’)dt’ 


(C7) 
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Substituting equations (C3), (C5), and (C6) into equation (C7) yields 


K = <(Tv) ne f t _ t n e [ Z “ v < t - t ’>5 *3 dt ’ 


(C8) 


by a similar argument, the attenuation coefficient for the downward -moving flux is given 

t>y 


^d 


w z r 

i ° L 


e d^d' t')dt' = (dv) ne f n e [z + v(t-t');t'] dt' (C9) 

c/t-t , 


The downward flux can be written in terms of this attenuation coefficient as 


*d = g 6XP 


Z-Z . 
min 


H 


o J 


~*r 


(CIO) 


The total particle flux at the point Z at time t is 


<p \i + (p d = n ( Z ’ = “ “i v exp 
2 




H 


o J 


(e'^Ve'S (Cll) 


The neutral density can therefore be written 


n(Z, t) = - n n exp 
2 1 


( Z - Z min) 


H 


o J 




(C 1 2) 


The net flux, with upward moving flux considered positive, is given by 


*n = *u " *d = j V exp 


(Z-Z . ) 
v min 7 


H 


o -J 


<> - e'S 


(C13) 
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The neutral particle continuity equation is identical to equation (B8) in this model. The 
divergence term, however, is 


d(p n 1 

-V • cp — n v exp 

n dZ 2 0 


(Z -WTe'*’ 1 - e~^ d . d »u _■% d *d ,-»dl (C14) 

H o JL H o dZ dZ j 


Performing the Z -derivatives of equations (C8) and (C14) and (C9), one obtains 


^d 

dZ 


= <ov> 


ne 


f 

Jt-t, 


dn e ( Z d ; » ^ne / 
e d dt’ — n (Z 


dZ 


e'“max’ ^ 


(Cl 5) 


and 




dZ 


(ov) 


ne 


f 

Jt-t 

u 


dn (Z t’) 
e v u’ 7 

dZ 


dt + 


<ov) 


ne 


n (Z . ; t-t ) 
e v min’ u' 


(C16) 


One cannot proceed further with the derivation without assuming a specific functional form 

for the altitude variation of n . I 

e 

creases exponentially with height, 


for the altitude variation of n g . If it is assumed that the electron number density m 


"e* 2 "') - "e.maxW ex " 


< Z max- Z > 


H, 


(C17) 


where H is the scale height of the electron layer, one obtains 


dn n 

e _ e 

dZ H 


(Cl 8) 


Incorporating equation (Cl 8) into equations (Cl 5) and (C16) allows them to be written 


d ^d _ ^d <° :v >ne n e (Z max’ t - t d) 


dZ H. 


(C19) 
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d K _ K < gv >ne 

dZ H v 
e 


n e^ Z min ; t_ V 


(C20) 


Substituting equations (C12), (C19), and (C20) into equation (C14) allows the diver- 
gence term to be written, 


-V. 0 = 


n(Z,t Xorv) 


ne 


(e'Se'S 


v(e^"-e'S v( V TU -^ d e' a ) 

— £ + + e n 




H (ov) no 
o ne 


H e<°^>ne 


e^min^-V 


+ e n e^ Z max ;t "V 


(C21) 


By hypothesis, the attenuation of the neutrals was very small for a traversal of the dis- 
charge, hence « 1, « 1. It will be shown later that for coherent oscillations 

to occur, one must have H g ~ H Q in the coherently oscillating region. Under these 
conditions, the first two terms within the brackets may be written 


Ai „“^d x v^ u e ^ u -i 


v(e -e u ) 
H o< OT >ne 


^o e 




(C22) 


H (crv) 
e ne 


2H ( av) 
o' ne 


If n g is an approximate spatial average of the electron number density, and T the tran- 
sit time across the oscillating region, one can write (crv) ne T > s0 that 

equation (C22) gives the approximate magnitude 


rfxjA- ib 2 ) n (Z -Z . )(i b,-\p ) 

yd ^u / ~ e v max mm AV a Y xx 


2H o<‘ ,v> ne 


2H. 


(C23) 


where one has made use of the relation 


Z -Z . = vT 
max mm 


(C24) 


Equation (C23) implies that, so long as the thickness of the oscillating layer 

Z -Z . is not very much larger than the scale height H , the first two terms in the 
max mm o 7 

brackets of equation (C21) are smaller than the last two by a factor of magnitude ^ 0 ~^ u - 
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Since this term is much less than unity by hypothesis, the first two terms will therefore 
be dropped from subsequent developments. 

Equation (C21) can then be written 


-V • cp n » n(Z,t)<CTv) ne 


r^^min^-y n e^ Z max’ t_ y 


1 + e 




1 + e 


(C25) 


Since « 1, one may expand the denominator to obtain 


- V * ^ \ n(Z ’ t)<crv> ne [“earnin’ Ut u} + n e (Z max^ 


(C26) 


One can now invoke the argument used previously to replace the time -retarded elec- 
tron number density with the time -averaged electron number density, since the velocity 
spread that is always present in actual cases will perform an averaging function. Experi- 
mental data relating to this assumption are presented in appendix E. Using a tilda to in- 
dicate a time average, we may write 


“earnin’ "" n e^ Z min^ 

(C27) 

n e^ Z max ; t-t d) "* n e^ Z max^ 

(C28) 


Equation (C26) can therefore be written 


-V . <p n « i n(Z, t)<ow) ne [n e (Z min ) + n e (Z max )] 


(C29) 


If we define a time and space averaged electron number density 


fi e = ^Pe (Z min) + %( Z max)] 


(C30) 


equation (C 29) can be written in terms of this average as 


-V- <p n = n(Z,t)n e <ov> ne 


(C31) 
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so that the neutral particle continuity equation may be written 

n= n(Z,t){(TV> ne [n e -n e (Z,t)] (C32) 

One can show by a physical argument virtually identical to that above that the electron 
continuity equation 


“e = ' V ‘ *ne + nn e <<rv> ne < C33 ) 

can be written in the symmetric form 

n e - “ n e (Z, t)(ow) ng [n-n(Z, t)^j ( c34 ) 

where n is the time- and space -averaged neutral number density. This physical argu- 
ment proceeds by replacing n g — n, and n — n e , and remembering that the signs of the 
arguments of the exponentials in equations (C4) and (CIO) are positive, since the ioniza- 
tion process enhances the electron density in the oscillating region. 

By following a derivation virtually identical to that in reference 11, one can show 
that when the peak-to-peak amplitude of n g (t) is small, the frequency of oscillation given 
by equations (C32) and (C34) is 


y'n(Z)n e (Z) (ov) ne 
277 


(C35) 


where n(Z) is the time -averaged neutral number density at the position Z, and n g (Z) 
in the time -averaged electron number density at the same position. In order for the gas 
to oscillate coherently and with the same frequency, it is necessary that the product of 
densities be independent of position, that is, 


n(Z)n e (Z) ~ constant (C36) 

Since the neutral density of planetary and stellar atmospheres falls off exponentially with 
height above the surface, equation (C36) implies that the electron number density must 
increase exponentially with altitude, and with the same scale height, over a significantly 
thick layer of the planetary or stellar atmosphere. Small departures from a constant 
product of n e n may be allowable if (orv) varies with height in such a way that the fre- 
quency (equation (C35)) is maintained constant. 
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APPENDIX D 


RELATION OF EXACT TO SMALL -AMPLITUDE SOLUTIONS 

The exact solutions discussed above are in many ways not as satisfactory as the 
small -amplitude solutions presented in references 10 , 11 , and 13 , since the latter gives 
the frequency G(ri) and the peak -to -peak amplitude as functions of 77 in closed form. 

When n eo is taken as the minimum value of n g (t), the parameter 77 used in the 
small amplitude solution is greater than unity, and the approximate solutions presented in 
reference 11 give a good qualitative account of the geometry and frequency of the wave- 
form of n g (t) . The peak -to -peak amplitude and G{ri) are given to within about 20 percent 
of the exact value by the small amplitude approximation. 

The agreement between the small -amplitude approximation and the exact solutions is 
much less satisfactory for 77 < 1 . 0 . If we denote the amplitude index used in the small- 

amplitude approximation by ri a , one can plot G(n_) and the peak-to-peak amplitude as a 

s s 

function of 77 for the relevant expressions in reference 11 . This has been done in fig- 
s 

ure 8 for the parameter G (77,,) and in figure 6 for the peak-to-peak amplitude. There is 

s 

fair agreement of the exact with the small -amplitude solution between 0. 80 ^ 77 ^ 1. 0, 

s 

but the two diverge significantly below j] a = 0 . 80 . This occurs because a mathematical 

s 

artifact of the small -amplitude approximation restricts the small -amplitude solutions to 

the range 2/3 < n ^ 1 . 0 , while the exact solutions range over 0 < < 1.0, 

The agreement of the exact with the small-amplitude solutions for 77 < 1 . 0 can be 

greatly improved by using the mapping between 77^ and 77. given on figure 21 . This 

mapping was obtained by requiring that the peak-to-peak amplitude of y(t) for the value 

of 77 e on the ordinate be equal to the peak-to-peak amplitude of y(t) obtained from the 

small -amplitude approximation for the value of 77. given on the abscissa. 

s 
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Electron amplitude index, r\ e 



Amplitude index, tj s - 2/3 


Figure 21. - Agreement between exact and small-amplitude solutions to equations (22) and (23) may be greatly 
improved by using heuristic transformation shewn, which assumes that exact and small-amplitude solutions 
have same peak-peak amplitude. For a particular value of electron amplitude index used in exact solution one 
should use equivalent value of amplitude index used in small-amplitude approximation. 
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APPENDIX E 


COHERENT OSCILLATION LIMITS AND NEUTRAL TRANSIT TIMES 

In the derivation in appendixes B and C, it is shown that the transit time delays, 
which occur as a neutral particle moves with finite velocity across the discharge, play a 
fundamental role in the oscillation mechanism. These time lags appeared in equations 
(B18) and (C26) of these appendixes, in which the divergence term is written 

-V • <p n « n(t)<av> ne n e (t-t Q ) (El) 

where t Q is the transit time delay associated with motion of a neutral particle of veloc- 
ity v across a characteristic discharge dimension R. One may define a phase angle 
/3 by nondimensionalizing the time t with respect to the period of oscillation T, 

- n e<^o> (E2) 

where 

£ sM= 2771/t (E 3 ) 

T 

and 

2nt rj 

R S ® = 2 TTV “ (E 4 ) 

° T v 

In appendixes B and C, it is argued that in actual applications, the velocity v is not 
a single value, but instead has a broad distribution of values. This in turn was used to 
justify replacing n e (t-t Q ) by the time -averaged electron density, n g . It is understandable 
that such an averaging would occur if ft s=a 27 r, since the transit time would then be com - 
parable to the period of oscillation. One might also expect the discharge to no longer 
oscillate coherently and with the same phase if /3 Q » 2tt, since the neutrals would then, 
on the average, take many periods to move across the discharge. One might then expect 
a wave of ionization to propagate across the discharge, much like a moving striation. The 
oscillations might also be expected to disappear if « 2ir, since the maximum possible 
time delays would then not be sufficient to perform an averaging function, and a detailed 
balance would exist between the divergence and loss terms of the continuity equation. 
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It is therefore of interest to determine the range of /3 over which coherent, in- 
phase oscillations can be expected to occur. The data reported in reference 13 was ac- 
cordingly reduced to determine the range of /3 for which coherent oscillations were 
observed in this experiment. The discharge radius was taken as R = 5. 0 centimeters, 
and the average velocity v also appearing in equation (E4) was taken to be the most pro- 
bable velocity of the molecules at a laboratory temperature of 283 K. The lowest and 
highest frequencies at which coherent in-phase oscillations were observed in deuterium, 
neon, and helium gas are listed in table in. 


TABLE III. - OBSERVED RANGE OF TRANSIT- TIME PHASE 
ANGLE FOR EXPERIMENTAL DATA REPORTED 
IN REFERENCE 15 


Species 

Lowest 
observed 
frequency . 
Hz 

Highest 

observed 

frequency. 

Hz 

Neutral 
transit 
time. 
;i sec 

Lower limit 
of phase 

,. ol 

Upper 
limit of 
phase 
a nide. 

,J oh- 

rad 

rad 

dc^ 

°2 

340 

50 500 

37. 7 j 

0. 0128 

4.6 

o 

o 

He 

106 

74 000 

84. 8 

.0090 

3.2 

6. 3 

He 

87 

22 500 

37. 7 

.0033 

1.2 

. 85 


Coherent oscillations may well have existed below the lowest frequencies shown, 
since these lower limits were determined by the disappearance of the light fluctuations 
into the background noise, and not by any loss of coherence. The upper limits were de- 
termined by the point above which propagating waves were observed in the discharge. 

The transit time of the neutrals across the discharge, in microseconds, is also given. 
The low and high values of P Q /2n are also indicated. 

It is, perhaps, surprising to find that the averaging process is effective when the 
phase angle /3 q j is as small as 1. 2°. Evidently therefore, only very small time lags 
are necessary to destroy the detailed balance between production and loss processes, and 
thereby permit the oscillation to take place. 
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